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LONG-TERM  GOAL 

The  ultimate  goal  of  this  research  is  the  development  of  a  complete  ice  dynamics  model  that  will 
include  lead  direction  and  ice  thickness  distribution  in  refrozen  leads. 

OBJECTIVES 

Northwest  Research  Associates,  Inc.  (NWRA)  had  two  objectives  during  this  past  year.  The  first  was 
to  use  stress  data  to  develop  a  new  formulation  of  ice  strength  to  be  used  in  the  anisotropic  model 
(Coon  et  al.,  1998).  The  new  formulation  was  developed  and  has  been  provided  to  Dr.  R.  S.  Pritchard 
for  his  use  in  model  development  for  PIPS3.0.  The  second  objective  for  this  year  was  to  review  other 
work  being  done  on  sea  ice  strength  in  Europe. 

APPROACH 

The  approach  this  year  was  to  use  measured  sea  ice  stress  and  models  to  formulate  the  behavior  of 
refrozen  leads.  Some  of  the  data  in  the  models  was  acquired  by  visits  to  several  European  universities. 
These  visits  were  accomplished  through  an  appointment  as  Adjunct  Scientist  to  the  International  ONR 
Office  in  London. 

WORK  COMPLETED 

The  work  with  SIMI  stress  data  analysis  was  completed  and  the  results  were  presented  in  a  paper, 
“Determination  of  Pack  Ice  Stress  from  Flatjack  Sensors,”  submitted  to  the  Journal  of  Cold  Regions 
Science  and  Technology.  Douglas  C.  Echert,  Gerald  S.  Knoke  and  myself  authored  this  paper. 
Important  insight  into  sea  ice  strength  was  gained  through  visits  to  the  Technical  Universities  of 
Denmark,  Finland,  and  Norway,  as  well  as  a  visit  to  the  Scott  Polar  Institute  of  Cambridge  University 
in  England.  On  each  visit,  I  made  a  presentation  on  the  development  of  the  Anisotropic  Sea  Ice  Model 
and  received  useful  review  of  the  model.  The  findings  from  the  visits  are  presented  in  the  section 
below.  The  results  of  this  year’s  work  were  also  presented  at  the  PIPS3.0  development  meeting  at  NIC 
in  July. 

RESULTS 

During  my  six-month  visit  to  the  Technical  University  of  Denmark,  I  worked  with  Leif  Toudal  on  a 
frazil/pancake  ice  model  for  the  Greenland  Sea  that  may  have  applications  in  support  of  PIPS3.0  for 
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the  marginal  ice  zones.  Other  European  travel  was  funded  through  the  International  ONR  Office  in 
London.  At  the  Technical  University  of  Finland,  they  are  working  on  ridge-building  forces  and  the 
strength  of  partly  consolidated  ridges.  The  work  on  ridge-building  forces  is  reported  by  Lensu  et  al. 
(1998)  and  may  be  particularly  important  for  ridges  in  the  marginal  ice  zone.  There  is  also  a  large- 
scale  field  experiment  determining  the  strength  of  first-year  ridges  by  utilizing  a  large-scale,  one  meter, 
downward  punch  test  through  the  ridge  keel.  The  results  may  be  of  value  in  determining  ridge  strength 
for  the  anisotropic  model  (Coon  et  al.,  1998).  There  is  work  at  the  University  of  Norway  on  the 
development  of  a  thermodynamic  model  for  the  consolidation  of  first-year  ridges  (see  Lpset  et  al., 
1998).  This  model  development  is  also  supported  by  annual  field  experiments.  The  results  from  this 
work  may  be  important  for  the  thermodynamic  consolidation  of  ridges  as  they  fit  into  an  anisotropic 
model.  At  the  Scott  Polar  Institute,  they  are  involved  in  a  European  project  called  Local  Ice  Cover 
Deformation  and  Mesoscale  Ice  Dynamics-Ice  State.  The  objectives  of  this  project  are: 

1.  To  describe  and  model  the  processes  involved  in  local  ice  cover  deformation,  such  as  ridging  and 
rafting,  and  to  incorporate  these  into  governing  equations  of  mesoscale  ice  dynamics;  and 

2.  To  identify  a  set  of  ice  cover  parameters  that  is  adequate  to  describe  an  ice  state,  which  can  ideally 
be  derived  from  remote  sensing  data  and  used  in  practical  navigation. 

The  results  from  this  project  may  aid  in  data  assimilation  for  the  PIPS3.0  model. 

For  the  anisotropic  model  of  Coon  et  al.  (1998),  it  is  important  to  develop  a  simple  relationship 
between  the  force  to  build  a  ridge  and  the  thickness  of  the  thinnest  ice  in  the  lead.  Considerable  model 
calculations  and  data  collection  related  to  this  problem  have  been  performed.  The  results  of  these 
investigations  can  be  summarized  in  a  simple  equation: 

F  =  110/r372  (  1  ) 

where  F  is  the  ridge-building  force  in  KN/m  and  h  is  the  ice  thickness  in  meters.  In  an  isotropic 
model,  F  would  be  the  lead  strength  and  h  the  thinnest  ice  in  the  lead.  This  equation  is  a  good 
approximation  to  the  average  ridge  building  force  and  the  viscoelastic -buckling  load  for  a  lead  (for 
loads  of  six-hour  duration).  This  is  in  good  agreement  with  Parmerter  and  Coon  (1973)  (for  friction  of 
0.3);  Coon  and  Lau  (1990);  Coon  et  al.  (1989a,  1989b);  and  Hopkins  (1998). 

The  drift  predicted  by  a  sea  ice  model  for  weak  ice  is  essentially  that  of  free  drift  where  ice  strength  is 
zero.  That  is  to  say,  the  strength  term  in  the  momentum  equation  is  small  by  comparison  to  other 
terms.  This  will  be  the  situation  as  the  ice  cover  forms  in  the  marginal  ice  zones.  Therefore,  it  seems 
reasonable  to  ask  the  question,  how  long  can  the  new  ice  in  the  marginal  ice  zone  be  modeled 
adequately  with  a  free-drift  dynamics  model?  One  can  get  an  estimate  of  this  time  by  examining  the 
growth  of  sea  ice  from  zero  to  some  thickness,  which  provides  a  ridging  strength  [given  in  Eq.  (1)] 
large  enough  to  be  important  in  the  sea  ice  dynamics  models. 

From  isotropic  ice  dynamics  model  calculations,  it  is  known  that  for  ice  strength  of  10  KN/m  the  ice 
motion  is  free  drift,  100  KN/m  is  strong  isotropic  ice  and  the  motion  departs  from  that  of  free  drift 
considerably.  In  an  anisotropic  model  there  are  two  strengths:  one  for  leads,  which  is  from  zero  to  the 
value  given  in  Eq.  (1)  and  a  second  strength  for  the  ice  that  contains  the  leads,  which  is  250  KN/m 
(Coon  et  al.,  1998).  It  would  appear  that  for  the  marginal  ice  zone  that  a  ridge  strength  between  Vi  and 


1  times  the  isotropic  strength  could  be  modeled  as  free  drift.  The  ice  thicknesses  associated  with  these 
strengths  are: 


Free  drift:  F  =  10  KN/m  and  Eq.  (1)  gives  h  =  .2  m 

Vi  of  isotropic  strength:  F  =  50  KN/m  and  Eq.  (1)  gives  h  =  .59  m 

Full  isotropic  strength:  F  =  100  KN/m  and  Eq.  (1)  gives  h  =  .94  m 

Attention  can  now  be  turned  to  the  time  period  required  to  grow  ice  to  thickness  of  .2  m,  .59  m,  and 
.94  m  in  the  marginal  ice  zone.  The  first  ice  to  form  is  frazil  and  pancake  ice.  Its  thickness  is 
determined  primarily  by  wind  speed  [see  Martin  and  Kauffman  (1981),  Bauer  and  Martin  (1983),  and 
Alam  and  Curry  (1998)].  In  the  work  on  frazil  and  pancake  ice  done  by  Toudal  and  Coon  (1999),  a 
simple  approximation  to  the  ice  thickness  can  be  given  by 

h  =  -\2  +  2.1Vci  (2) 

where  Va  is  the  10  m  wind  speed  in  m/sec  for  Va  >4.3. 

After  an  ice  cover  has  formed,  the  ice  growth  is  governed  by  cumulative  freezing  days  in  °C  as  given 
by  Zubov  (1944)  to  be: 


h2  +  50 h  -  8R  =  0  (3) 

where  R  is  cumulative  freezing  days  in  °C. 

As  an  example,  it  will  be  assumed  that  when  the  ice  cover  forms,  there  is  Va  =  12  and  thereafter,  the 
temperature  is  -20°C.  Then  the  frazil  and  pancake  ice  would  first  form  to  a  thickness  of  0.2  m  as  given 
by  Eq.  (2).  The  question  now  becomes,  how  long  will  it  take  to  grow  the  ice  to  .59  m  and  .94  m. 

From  Eq.  (3)(with  a  temperature  of  -20°C)  it  will  take: 

9  days  to  grow  .2  m  of  ice; 

40  days  to  grow  .59  m  of  ice;  and 
84  days  to  grow  .94  m  of  ice. 

It  can  now  be  seen  that  in  this  example,  the  .59  m  of  ice  can  be  grown  in  40-9  =  31  days  and  the  .94  m 
of  ice  can  be  grown  in  84-9  =  75  days.  The  first  .2  m  is  grown  by  the  action  of  the  wind  [Eq.  (2)]  so 
the  nine  days  required  to  grow  it  by  freezing  can  be  subtracted  from  the  total  number  of  days  that 
would  be  required  to  form  it  from  freezing. 

The  conclusion  from  this  is:  The  ice  in  the  marginal  ice  zone  can  be  modeled  as  free  drift  for  the  first 
one  to  two  months. 
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